Multiplicity of Ordered Phases in Frustrated Systems 
Obtained from Hard-Spin Mean-Field Theory 



O 

o 
o 

(N 

» 

o 
O 

(N 



CO 



I 

C 

O 

o 



> 

00 



o 
o 



c 

o 
o 



X 
S3 



Hiiseyin Kaya 1,2 and A. Nihat Berker 1>2 ' 3 
Feza Giirsey Research Institute for Basic Sciences, 81220 Cengelkdy, Istanbul, Turkey 
2 Department of Physics, Istanbul Technical University, 80626 Maslak, Istanbul, Turkey 
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 

Random quenched dilution of the triangular-lattice antiferromagnetic Ising model locally relieves 
frustration, leading to ordering phenomena. We have studied this system, under such dilution of 
one sublattice, using hard-spin mean-field theory. After a threshold dilution, two sublattices develop 
non-zero magnetizations of equal magnitude and opposite signs, as all three sublattices exhibit spin- 
glass order. In this phase, multiple sets of ordered solutions occur. A phase diagram is obtained in 
dilution fraction and temperature. 
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A " rugged free-energy landscape" is often mentioned as 
a distinctive characteristic in the discussions of spin-glass 
systems jjj. Concrete support for such a phenomenon 
would derive from multiple solutions, not related to each 
other by a global symmetry, of self-consistent order- 
parameter equations. This has not been previously ob- 
tained for any system with a realistic spatial connectivity. 
In this work, we do find such multiple solutions, not re- 
lated by global symmetry, in a random frustrated system 
with realistic spatial connectivity, namely the quench- 
diluted triangular-lattice antiferromagnetic Ising model, 
studied via the closed-form implementation of hard-spin 
mean- field theory @-|5J. 

The antiferromagnetic Ising model, with Hamiltonian 
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J > 0, 



(1) 



where Sj = ±1 at each site i of a triangular lattice and 
< ij > indicates summation over nearest-neighbor pairs 
of sites, is fully frustrated [13]: In each elementary trian- 
gle, one of the three nearest-neighbor antiferromagnetic 
interactions is dissatisfied when the energy is minimized. 
This leads, macroscopically, to a highly degenerate sys- 
tem that is disordered at all non-zero (1/J > 0) tem- 
peratures Hj. Random quenched dilution of the sys- 
tem relieves the frustration at random localities and can 
be expected to lead to ordering phenomena. Indeed, a 
Monte Carlo study with random quenched dilution of all 
three sublattices equivalently has indicated spin-glass or- 
der pq| . We consider the random quenched dilution of 
the sites of one of three sublattices. After a threshold 
dilution, the system exhibits uniform and opposite mag- 
netizations in two sublattices and spin-glass order, i.e., 
spins frozen in random directions, in the quench-dilutcd 
sublattice. A phase diagram is obtained in the variables 
of dilution fraction and temperature. Within the ordered 
phase, for a fixed dilution fraction and temperature, a 
multiplicity of solutions, distinguished by different val- 
ues of the local and global order parameters, is obtained 
to the hard-spin mean- field equations. 



We use hard-spin mean-field theory, a method which is 
almost as simply enunciated as usual mean-field theory, 
but which conserves frustration [p|-[l2||. Consequently, 
this method has yielded, for example, the lack of or- 
der in the undiluted triangular-lattice antiferromagnetic 
Ising model and the ordering that occurs when a uniform 
magnetic field is applied to the system, in a quantita- 
tively accurate phase diagram in the temperature versus 
magnetic field variables Hard-spin mean-field 

theory also yields the lack of order in the one-dimensional 
Ising ferromagnet and the occurrence of order in the two- 
dimensional Ising ferromagnet, the latter with an onset 
temperature improved over usual mean-field theory [0. 
Hard-spin mean-field theory has also been successfully 
applied to complicated systems that exhibit a variety 
of ordering behaviors, such as three-dimensional stacked 
frustrated systems |2||| and higher-spin systems [Q. 

The self-consistent equation for local magnetizations 
in hard-spin mean-field theory is 



= Y1 Y[p( m f' s j) tanh(-J^\ 

{«} j J j' 



(2) 



where the product over j and sum over j' run over all 
non-diluted sites neighboring site i, and the single-site 
spin probability distribution p(rnj\Sj) is (1 + rrijSj)/2. 
The outer summation is over the ±1 values of the spins 
at the undiluted sites neighboring site i. Thus, the spin 
at each site feels the anti-aligning field due to the full 
(i.e., hard-) spin of each of its neighbors. Eq.(^) is a set 
of coupled equations for all the local magnetizations and 
is solved iteratively for a given realization of dilution in 
a finite but large system. 

Alternatively, a further approximation is to impose 
sublatticewise uniformity, rrij — m a for each sublattice 
a = a, b, c, and to average the self-consistent equation 
over all realizations of quenched site dilution, 
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In Eq.(j^), the parentheses enclose the right-hand side of 
Eq.(|J). This is summed over the 2 6 possible quenched 
environments {77} of site i. Each quenched environment 
{77} occurs with a probability Q composed of six factors 
qj, with qj = (1 — pj) for each quench-diluted neighbor 
j and qj — pj for each undiluted neighbor j. Eq.(||) can 
be compactly rewritten as 



j 773=0,1 



=±1 



(4) 

where the product over j and sum over j' run over all 
sites neighboring site i, the single-site quenched-dilution 
probability distribution q(pj] r/j) is 1 — pj — r/j(l — 2pj), 
and Pj = p a for each sublattice a = a,b,c. Eq.(||) or, 
equivalently, Eq.(|j) is solved for (m a , mi,, m c ) for given 

(Pa,Pb,Pc) ■ 

Whereas hard-spin mean-field theory [Eq.(|])] yields 
the variations in the local magnetizations within each 
sublattice due to differently quenched local environ- 
ments, the implementation of Eq.(^) is a further ap- 
proximation over hard-spin mean-field theory: While still 
incorporating frustration, it imposes sublatticewise uni- 
form magnetizations. Eq.Q is a set of three coupled 
equations, whereas Eq.(||) is a set of N coupled equa- 
tions, where N is the size of the system. 

Upon random quenched dilution of one sublattice, 
the frustrated triangular-lattice Ising model does in- 
deed show long-range order, as for example depicted in 
Figs.l. The two undiluted sublattices (labeled a and 6), 
which are now subject to random unfrustrated localities 
at the dilution points of the other sublattice (labeled 
c), develop non-zero sublattice-averaged magnetizations, 
m a — —rut, at low temperatures. For low dilutions, these 
magnetizations [Fig. 1(a)] show an initial slow growth at 
onset as temperature is lowered and do not saturate at 
zero temperature. The diluted sublattice c is also sub- 
ject to local liftings of frustration, due to the spatially 
non-uniform magnetizations of sublattices a and b, but 
this is a secondary and, therefore, weaker effect, and sub- 
lattice c does not develop a non-zero sublattice-averaged 
magnetization, m c = 0. 

All three sublattices develop, within the ordered phase, 
non-zero spin-glass order , i.e. , randomly frozen order, 
with order parameters 
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where N a is the number of spins of sublattice a. 
that, for the quenched-diluted sublattice, 



(5) 



Note 
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At high dilutions, the spin- glass ordering trend shows 
reentrance (as temperature is lowered, increases and 
then decreases) on the undiluted sublattices [Fig. 1(b)] 
and double reentrance (increase, decrease, and again in- 
crease) on the diluted sublattice [Fig. 1(c)]. Maximal 
zero-temperature spin-glass order occurs at intermediate 
dilutions, as seen in Figs.2(b,c). 
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FIG. 1. Finite-temperature order in the random quench- 
diluted triangular-lattice antiferromagnetic Ising model: 
(a) Magnetizations of the two undiluted sublattices from 
hard-spin mean-field theory (full curves). The three curves 
correspond to p = 0.75, 0.5, 0.09375, which can be iden- 
tified by the respectively increasing critical temperatures, 
i.e., increasing x-axis intercepts. The result from the fur- 
ther approximation of Eq.(3) is given with the dashed 
curve. The quench-diluted sublattice has zero magne- 
tization, (b) Spin-glass order parameter of the undi- 
luted sublattices. The curves, again distinguished by 
their respectively increasing critical temperatures, are for 
p = 0.9375, 0.890625, 0.75, 0.625, 0.5, 0.140625, 0.09375. Note 
the reentrant behavior in the magnitudes, for high dilutions, 
(c) Spin-glass order parameter of the quench-diluted sublat- 
tice. The values of p are as in (b). Here, the spin-glass order 
for p = 0.09375 is away from zero only at higher tempera- 
tures. The spin-glass order for p = 0.140625 exhibits doubly 
reentrant behavior. The full curves in Fig.l are for a fixed 
realization of the quenched disorder in a 24 x 24 system. 



2 




O 0.5 1 

Sublattice occupation, p c 

FIG. 2. Zero-temperature sublattice magnetizations and 
spin-glass order parameters from hard-spin mean-field the- 
ory. The sublattice magnetizations (a) do not saturate at 
zero-temperature, except for the full dilution (hexagonal lat- 
tice) limit. The dashed curve in (a) is the result of the further 
approximation of Eq.(3). The zero-temperature spin-glass or- 
der (b,c) is maximal at intermediate dilutions. Multiple solu- 
tions of the hard-spin mean-field theory equations are seen, in 
addition to the most stable set in depicted in Figs.l. (For each 
solution, the same symbol is used in Figs. 2-6). This figure is 
obtained for 1/J = 0.0001, in a 30 x 30 system. 

The phase diagram of the system is shown in Fig. 3. It 
is seen that a threshold dilution of 0.042 of one sublattice 
is needed for ordering, i.e., the occupancy p has to be be- 
low 0.958 for ordering. Fig. 3 shows the phase boundaries 
obtained for two realizations of quenched dilutions of a 
24 x 24 system and the result from averaging over 15 such 
realizations. Also shown in dashed in Fig. 3 is the further 



approximation of Eq.< 
obeys the equation 



i). This dashed phase boundary 
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where t n = tanh(nJ), 
and gives a dilution threshold of 0.125. 
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FIG. 3. Phase diagram of the random quench-diluted tri- 
angular-lattice antiferromagnetic Ising model. The full curves 
show the phase boundaries obtained by hard-spin mean-field 
theory for two realizations of quenched dilutions of a 24 x 24 
system. The losanges show the result from averaging over 15 
such realizations. The dashed curve is the result from the fur- 
ther approximation of Eq.(3); this curve is given analytically 
by Eq.(6). 



3p(l-p) 2 /i + (l-p) 3 /o = 4/3 (7) 
6t 4 + 5t 2 )/8, 
3t 3 + 2*0/4, 
2t 2 )/2, 




FIG. 4. 
= 0.75. 
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Finite-temperature multiplicity of solutions for 
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Within the ordered phase, we find that the hard-spin 
mean-field equations (||) admit, as seen in Figs. 2, 4, a mul- 
tiplicity of solutions, in addition to the set depicted in 
Figs.l. The latter is the most stable solution, in the 
sense that it has the largest basin of attraction under the 
iterative solution of Eq. ^ . The other solutions appear 
at different temperatures below the onset temperature 
for the most stable solution. Since, in this study, the 
number of the sets of solutions increased in going from 
the 24 x 24 system to the 30 x 30 system (depicted in 
Figs. 2, 4), it can be inferred that the solutions become 
numerous in the infinite system. 



Figs. 5, 6 depict the overlaps between pairs of solutions 
(A,B), 
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FIG. 5. Overlap Q between the different solutions [see 
Eq.(7a)] on the undiluted sublattices, for p = 0.75. The arrow 
indicates the onset of order. 
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FIG. 6. Overlap A between the different solutions [see 
Eq.(7b)] on the undiluted sublattices, for p — 0.75 . The 
arrow indicates the onset of order. In inset schematizes the 
separations between the different solutions, showing no ultra- 
metricity. 
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and the distances 
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From A^ B in Fig. 6, we can deduce the separation, in 
local-order-parameter space, between the different solu- 
tions. This is shown schematically in the inset of the 
figure. We note that the different solutions are not ul- 
trametrically Jl7j related, since no isosceles-triangle rela- 
tions are seen. Thus, it may well be that ultrametricity 
is a property specific to the infinitely connected lattice. 
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